Tapcin, an In Vivo Active Dual Topoisomerase I/II Inhibitor Discovered by Synthetic Bioinformatic Natural Product (Syn-BNP)-Coupled Metagenomics.

DNA topoisomerases are attractive targets for anticancer agents. Dual topoisomerase I/II inhibitors are particularly appealing due to their reduced rates of resistance. A number of therapeutically relevant topoisomerase inhibitors are bacterial natural products. Mining the untapped chemical diversity encoded by soil microbiomes presents an opportunity to identify additional natural topoisomerase inhibitors. Here we couple metagenome mining, bioinformatic structure prediction algorithms, and chemical synthesis to produce the dual topoisomerase inhibitor tapcin. Tapcin is a mixed p-aminobenzoic acid (PABA)-thiazole with a rare tri-thiazole substructure and picomolar antiproliferative activity. Tapcin reduced colorectal adenocarcinoma HT-29 cell proliferation and tumor volume in mouse hollow fiber and xenograft models, respectively. In both studies it showed similar activity to the clinically used topoisomerase I inhibitor irinotecan. The study suggests that the interrogation of soil microbiomes using synthetic bioinformatic natural product methods has the potential to be a rewarding strategy for identifying potent, biomedically relevant, antiproliferative agents.

High-throughput retrieval of target sequences from complex clone libraries using CRISPRi.

The capture of metagenomic DNA in large clone libraries provides the opportunity to study microbial diversity that is inaccessible using culture-dependent methods. In this study, we harnessed nuclease-deficient Cas9 to establish a CRISPR counter-selection interruption circuit (CCIC) that can be used to retrieve target clones from complex libraries. Combining modern sequencing methods with CCIC cloning allows for rapid physical access to the genetic diversity present in natural ecosystems.

Discovery of Paenibacillaceae Family Gram-Negative-Active Cationic Lipopeptide Antibiotics Using Evolution-Guided Chemical Synthesis.

Cationic nonribosomal lipopeptides (CNRLPs) from Paenibacillus spp. have been a rewarding source of Gram-negative-active antibiotics. Here we systematically screened sequenced bacterial genomes for CNRLP biosynthetic gene clusters (BGCs) that we predicted might encode additional Gram-negative-active antibiotics. Total chemical synthesis of the bioinformatically predicted products of seven such BGCs led to our identification of new laterocidine, tridecaptin, and paenibacterin-like antibiotics with potent activity against both multiple-drug-resistant Gram-negative and Gram-positive pathogens.

Bioinformatic prospecting and synthesis of a bifunctional lipopeptide antibiotic that evades resistance.

Emerging resistance to currently used antibiotics is a global public health crisis. Because most of the biosynthetic capacity within the bacterial kingdom has remained silent in previous antibiotic discovery efforts, uncharacterized biosynthetic gene clusters found in bacterial genome-sequencing studies remain an appealing source of antibiotics with distinctive modes of action. Here, we report the discovery of a naturally inspired lipopeptide antibiotic called cilagicin, which we chemically synthesized on the basis of a detailed bioinformatic analysis of the cil biosynthetic gene cluster. Cilagicin's ability to sequester two distinct, indispensable undecaprenyl phosphates used in cell wall biosynthesis, together with the absence of detectable resistance in laboratory tests and among multidrug-resistant clinical isolates, makes it an appealing candidate for combating antibiotic-resistant pathogens.

Lapcin, a potent dual topoisomerase I/II inhibitor discovered by soil metagenome guided total chemical synthesis.

In natural product discovery programs, the power of synthetic chemistry is often leveraged for the total synthesis and diversification of characterized metabolites. The synthesis of structures that are bioinformatically predicted to arise from uncharacterized biosynthetic gene clusters (BGCs) provides a means for synthetic chemistry to enter this process at an early stage. The recent identification of non-ribosomal peptides (NRPs) containing multiple ρ-aminobenzoic acids (PABAs) led us to search soil metagenomes for BGCs that polymerize PABA. Here, we use PABA-specific adenylation-domain sequences to guide the cloning of the lap BGC directly from soil. This BGC was predicted to encode a unique N-acylated PABA and thiazole containing structure. Chemical synthesis of this structure gave lapcin, a dual topoisomerase I/II inhibitor with nM to pM IC50s against diverse cancer cell lines. The discovery of lapcin highlights the power of coupling metagenomics, bioinformatics and total chemical synthesis to unlock the biosynthetic potential contained in even complex uncharacterized BGCs.

A naturally inspired antibiotic to target multidrug-resistant pathogens.

Gram-negative bacteria are responsible for an increasing number of deaths caused by antibiotic-resistant infections1,2. The bacterial natural product colistin is considered the last line of defence against a number of Gram-negative pathogens. The recent global spread of the plasmid-borne mobilized colistin-resistance gene mcr-1 (phosphoethanolamine transferase) threatens the usefulness of colistin3. Bacteria-derived antibiotics often appear in nature as collections of similar structures that are encoded by evolutionarily related biosynthetic gene clusters. This structural diversity is, at least in part, expected to be a response to the development of natural resistance, which often mechanistically mimics clinical resistance. Here we propose that a solution to mcr-1-mediated resistance might have evolved among naturally occurring colistin congeners. Bioinformatic analysis of sequenced bacterial genomes identified a biosynthetic gene cluster that was predicted to encode a structurally divergent colistin congener. Chemical synthesis of this structure produced macolacin, which is active against Gram-negative pathogens expressing mcr-1 and intrinsically resistant pathogens with chromosomally encoded phosphoethanolamine transferase genes. These Gram-negative bacteria include extensively drug-resistant Acinetobacter baumannii and intrinsically colistin-resistant Neisseria gonorrhoeae, which, owing to a lack of effective treatment options, are considered among the highest level threat pathogens4. In a mouse neutropenic infection model, a biphenyl analogue of macolacin proved to be effective against extensively drug-resistant A. baumannii with colistin-resistance, thus providing a naturally inspired and easily produced therapeutic lead for overcoming colistin-resistant pathogens.

Identification of structurally diverse menaquinone-binding antibiotics with in vivo activity against multidrug-resistant pathogens.

The emergence of multidrug-resistant bacteria poses a threat to global health and necessitates the development of additional in vivo active antibiotics with diverse modes of action. Directly targeting menaquinone (MK), which plays an important role in bacterial electron transport, is an appealing, yet underexplored, mode of action due to a dearth of MK-binding molecules. Here we combine sequence-based metagenomic mining with a motif search of bioinformatically predicted natural product structures to identify six biosynthetic gene clusters that we predicted encode MK-binding antibiotics (MBAs). Their predicted products (MBA1-6) were rapidly accessed using a synthetic bioinformatic natural product approach, which relies on bioinformatic structure prediction followed by chemical synthesis. Among these six structurally diverse MBAs, four make up two new MBA structural families. The most potent member of each new family (MBA3, MBA6) proved effective at treating methicillin-resistant Staphylococcus aureus infection in a murine peritonitis-sepsis model. The only conserved feature present in all MBAs is the sequence 'GXLXXXW', which we propose represents a minimum MK-binding motif. Notably, we found that a subset of MBAs were active against Mycobacterium tuberculosis both in vitro and in macrophages. Our findings suggest that naturally occurring MBAs are a structurally diverse and untapped class of mechanistically interesting, in vivo active antibiotics.

Genome-wide characterization of human minisatellite VNTRs: population-specific alleles and gene expression differences.

Variable Number Tandem Repeats (VNTRs) are tandem repeat (TR) loci that vary in copy number across a population. Using our program, VNTRseek, we analyzed human whole genome sequencing datasets from 2770 individuals in order to detect minisatellite VNTRs, i.e., those with pattern sizes ≥7 bp. We detected 35 638 VNTR loci and classified 5676 as commonly polymorphic (i.e. with non-reference alleles occurring in >5% of the population). Commonly polymorphic VNTR loci were found to be enriched in genomic regions with regulatory function, i.e. transcription start sites and enhancers. Investigation of the commonly polymorphic VNTRs in the context of population ancestry revealed that 1096 loci contained population-specific alleles and that those could be used to classify individuals into super-populations with near-perfect accuracy. Search for quantitative trait loci (eQTLs), among the VNTRs proximal to genes, indicated that in 187 genes expression differences correlated with VNTR genotype. We validated our predictions in several ways, including experimentally, through the identification of predicted alleles in long reads, and by comparisons showing consistency between sequencing platforms. This study is the most comprehensive analysis of minisatellite VNTRs in the human population to date.

Biosynthetic Interrogation of Soil Metagenomes Reveals Metamarin, an Uncommon Cyclomarin Congener with Activity against Mycobacterium tuberculosis.

Tuberculosis (TB) remains one of the deadliest infectious diseases. Unfortunately, the development of antibiotic resistance threatens our current therapeutic arsenal, which has necessitated the discovery and development of novel antibiotics against drug-resistant Mycobacterium tuberculosis (Mtb). Cyclomarin A and rufomycin I are structurally related cyclic heptapeptides assembled by nonribosomal peptide synthetases (NRPSs), which show potent anti-Mtb activity with a new cellular target, the caseinolytic protein ClpC1. An NRPS adenylation domain survey using DNA extracted from ∼2000 ecologically diverse soils found low cyclomarin/rufomycin biosynthetic diversity. In this survey, a family of cyclomarin/rufomycin-like biosynthetic gene clusters (BGC) that encode metamarin, an uncommon cyclomarin congener with potent activity against both Mtb H37Rv and multidrug-resistant Mtb clinical isolates was identified. Metamarin effectively inhibits Mtb growth in murine macrophages and increases the activities of ClpC1 ATPase and the associated ClpC1/P1/P2 protease complex, thus causing cell death by uncontrolled protein degradation.

Sleep-wake disorders in Alzheimer's disease: further genetic analyses in relation to objective sleep measures.

This paper presents updated analyses on the genetic associations of sleep disruption in individuals with Alzheimer's disease (AD). We published previously a study of the association between single nucleotide polymorphisms (SNPs) found in eight genes related to circadian rhythms and objective measures of sleep-wake disturbances in 124 individuals with AD. Here, we present new relevant analyses using polygenic risk scores (PRS) and variable number tandem repeats (VNTRs) enumerations. PRS were calculated using the genetic data from the original participants and relevant genome wide association studies (GWAS). VNTRs for the same circadian rhythm genes studied with SNPs were obtained from a separate cohort of participants using whole genome sequencing (WGS). Objectively (wrist actigraphy) determined wake after sleep onset (WASO) was used as a measure of sleep disruption. None of the PRS were associated with sleep disturbance. Computer analyses using VNTRseek software generated a total of 30 VNTRs for the circadian-related genes but none appear relevant to our objective sleep measure. In addition, of 71 neurotransmitter function-related genes, 29 genes had VNTRs that differed from the reference VNTR, but it was not clear if any of these might affect circadian function in AD patients. Although we have not found in either the current analyses or in our previous published analyses of SNPs any direct linkages between identified genetic factors and WASO, research in this area remains in its infancy.

BpWrapper: BioPerl-based sequence and tree utilities for rapid prototyping of bioinformatics pipelines.

BACKGROUND: Automated bioinformatics workflows are more robust, easier to maintain, and results more reproducible when built with command-line utilities than with custom-coded scripts. Command-line utilities further benefit by relieving bioinformatics developers to learn the use of, or to interact directly with, biological software libraries. There is however a lack of command-line utilities that leverage popular Open Source biological software toolkits such as BioPerl ( http://bioperl.org ) to make many of the well-designed, robust, and routinely used biological classes available for a wider base of end users. RESULTS: Designed as standard utilities for UNIX-family operating systems, BpWrapper makes functionality of some of the most popular BioPerl modules readily accessible on the command line to novice as well as to experienced bioinformatics practitioners. The initial release of BpWrapper includes four utilities with concise command-line user interfaces, bioseq, bioaln, biotree, and biopop, specialized for manipulation of molecular sequences, sequence alignments, phylogenetic trees, and DNA polymorphisms, respectively. Over a hundred methods are currently available as command-line options and new methods are easily incorporated. Performance of BpWrapper utilities lags that of precompiled utilities while equivalent to that of other utilities based on BioPerl. BpWrapper has been tested on BioPerl Release 1.6, Perl versions 5.10.1 to 5.25.10, and operating systems including Apple macOS, Microsoft Windows, and GNU/Linux. Release code is available from the Comprehensive Perl Archive Network (CPAN) at https://metacpan.org/pod/Bio::BPWrapper . Source code is available on GitHub at https://github.com/bioperl/p5-bpwrapper . CONCLUSIONS: BpWrapper improves on existing sequence utilities by following the design principles of Unix text utilities such including a concise user interface, extensive command-line options, and standard input/output for serialized operations. Further, dozens of novel methods for manipulation of sequences, alignments, and phylogenetic trees, unavailable in existing utilities (e.g., EMBOSS, Newick Utilities, and FAST), are provided. Bioinformaticians should find BpWrapper useful for rapid prototyping of workflows on the command-line without creating custom scripts for comparative genomics and other bioinformatics applications.

BitPAl: a bit-parallel, general integer-scoring sequence alignment algorithm.

MOTIVATION: Mapping of high-throughput sequencing data and other bulk sequence comparison applications have motivated a search for high-efficiency sequence alignment algorithms. The bit-parallel approach represents individual cells in an alignment scoring matrix as bits in computer words and emulates the calculation of scores by a series of logic operations composed of AND, OR, XOR, complement, shift and addition. Bit-parallelism has been successfully applied to the longest common subsequence (LCS) and edit-distance problems, producing fast algorithms in practice. RESULTS: We have developed BitPAl, a bit-parallel algorithm for general, integer-scoring global alignment. Integer-scoring schemes assign integer weights for match, mismatch and insertion/deletion. The BitPAl method uses structural properties in the relationship between adjacent scores in the scoring matrix to construct classes of efficient algorithms, each designed for a particular set of weights. In timed tests, we show that BitPAl runs 7-25 times faster than a standard iterative algorithm. AVAILABILITY AND IMPLEMENTATION: Source code is freely available for download at http://lobstah.bu.edu/BitPAl/BitPAl.html. BitPAl is implemented in C and runs on all major operating systems. CONTACT: jloving@bu.edu or yhernand@bu.edu or gbenson@bu.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

VNTRseek-a computational tool to detect tandem repeat variants in high-throughput sequencing data.

DNA tandem repeats (TRs) are ubiquitous genomic features which consist of two or more adjacent copies of an underlying pattern sequence. The copies may be identical or approximate. Variable number of tandem repeats or VNTRs are polymorphic TR loci in which the number of pattern copies is variable. In this paper we describe VNTRseek, our software for discovery of minisatellite VNTRs (pattern size ≥ 7 nucleotides) using whole genome sequencing data. VNTRseek maps sequencing reads to a set of reference TRs and then identifies putative VNTRs based on a discrepancy between the copy number of a reference and its mapped reads. VNTRseek was used to analyze the Watson and Khoisan genomes (454 technology) and two 1000 Genomes family trios (Illumina). In the Watson genome, we identified 752 VNTRs with pattern sizes ranging from 7 to 84 nt. In the Khoisan genome, we identified 2572 VNTRs with pattern sizes ranging from 7 to 105 nt. In the trios, we identified between 2660 and 3822 VNTRs per individual and found nearly 100% consistency with Mendelian inheritance. VNTRseek is, to the best of our knowledge, the first software for genome-wide detection of minisatellite VNTRs. It is available at http://orca.bu.edu/vntrseek/.

Inter- and intra-specific pan-genomes of Borrelia burgdorferi sensu lato: genome stability and adaptive radiation.

BACKGROUND: Lyme disease is caused by spirochete bacteria from the Borrelia burgdorferi sensu lato (B. burgdorferi s.l.) species complex. To reconstruct the evolution of B. burgdorferi s.l. and identify the genomic basis of its human virulence, we compared the genomes of 23 B. burgdorferi s.l. isolates from Europe and the United States, including B. burgdorferi sensu stricto (B. burgdorferi s.s., 14 isolates), B. afzelii (2), B. garinii (2), B. "bavariensis" (1), B. spielmanii (1), B. valaisiana (1), B. bissettii (1), and B. "finlandensis" (1). RESULTS: Robust B. burgdorferi s.s. and B. burgdorferi s.l. phylogenies were obtained using genome-wide single-nucleotide polymorphisms, despite recombination. Phylogeny-based pan-genome analysis showed that the rate of gene acquisition was higher between species than within species, suggesting adaptive speciation. Strong positive natural selection drives the sequence evolution of lipoproteins, including chromosomally-encoded genes 0102 and 0404, cp26-encoded ospC and b08, and lp54-encoded dbpA, a07, a22, a33, a53, a65. Computer simulations predicted rapid adaptive radiation of genomic groups as population size increases. CONCLUSIONS: Intra- and inter-specific pan-genome sizes of B. burgdorferi s.l. expand linearly with phylogenetic diversity. Yet gene-acquisition rates in B. burgdorferi s.l. are among the lowest in bacterial pathogens, resulting in high genome stability and few lineage-specific genes. Genome adaptation of B. burgdorferi s.l. is driven predominantly by copy-number and sequence variations of lipoprotein genes. New genomic groups are likely to emerge if the current trend of B. burgdorferi s.l. population expansion continues.

Pervasive recombination and sympatric genome diversification driven by frequency-dependent selection in Borrelia burgdorferi, the Lyme disease bacterium.

How genomic diversity within bacterial populations originates and is maintained in the presence of frequent recombination is a central problem in understanding bacterial evolution. Natural populations of Borrelia burgdorferi, the bacterial agent of Lyme disease, consist of diverse genomic groups co-infecting single individual vertebrate hosts and tick vectors. To understand mechanisms of sympatric genome differentiation in B. burgdorferi, we sequenced and compared 23 genomes representing major genomic groups in North America and Europe. Linkage analysis of >13,500 single-nucleotide polymorphisms revealed pervasive horizontal DNA exchanges. Although three times more frequent than point mutation, recombination is localized and weakly affects genome-wide linkage disequilibrium. We show by computer simulations that, while enhancing population fitness, recombination constrains neutral and adaptive divergence among sympatric genomes through periodic selective sweeps. In contrast, simulations of frequency-dependent selection with recombination produced the observed pattern of a large number of sympatric genomic groups associated with major sequence variations at the selected locus. We conclude that negative frequency-dependent selection targeting a small number of surface-antigen loci (ospC in particular) sufficiently explains the maintenance of sympatric genome diversity in B. burgdorferi without adaptive divergence. We suggest that pervasive recombination makes it less likely for local B. burgdorferi genomic groups to achieve host specialization. B. burgdorferi genomic groups in the northeastern United States are thus best viewed as constituting a single bacterial species, whose generalist nature is a key to its rapid spread and human virulence.

Linkage disequilibrium based genotype calling from low-coverage shotgun sequencing reads.

BACKGROUND: Recent technology advances have enabled sequencing of individual genomes, promising to revolutionize biomedical research. However, deep sequencing remains more expensive than microarrays for performing whole-genome SNP genotyping. RESULTS: In this paper we introduce a new multi-locus statistical model and computationally efficient genotype calling algorithms that integrate shotgun sequencing data with linkage disequilibrium (LD) information extracted from reference population panels such as Hapmap or the 1000 genomes project. Experiments on publicly available 454, Illumina, and ABI SOLiD sequencing datasets suggest that integration of LD information results in genotype calling accuracy comparable to that of microarray platforms from sequencing data of low-coverage. A software package implementing our algorithm, released under the GNU General Public License, is available at http://dna.engr.uconn.edu/software/GeneSeq/. CONCLUSIONS: Integration of LD information leads to significant improvements in genotype calling accuracy compared to prior LD-oblivious methods, rendering low-coverage sequencing as a viable alternative to microarrays for conducting large-scale genome-wide association studies.

Fast, adaptive evolution at a bacterial host-resistance locus: the PFam54 gene array in Borrelia burgdorferi.

Microbial pathogens have evolved sophisticated mechanisms for evasion of host innate and adaptive immunities. PFam54 is the largest paralogous gene family in the genomes of Borrelia burgdorferi, the Lyme disease bacterium. One member of PFam54, the complement-regulator acquiring surface proteins 1 (BbCrasp-1), is able to abort the alternative pathway of complement activation via binding human complement-regulator factor H (FH). The gene coding for BbCRASP-1 exists in a tandem array of PFam54 genes in the B. burgdorferi genome, a result apparently of repeated gene duplications. To help elucidate the functions of the large number of PFam54 genes, we performed phylogenomic and structural analyses of the PFam54 gene array from ten B. burgdorferi genomes. Analyses based on gene tree, genome synteny, and structural models revealed rapid adaptive evolution of this array through gene duplication, gene loss, and functional diversification. Individual PFam54 genes, however, do not show high intra-population sequence polymorphisms as genes providing evasion from adaptive immunity generally do. PFam54 members able to bind human FH are not monophyletic, suggesting that human FH affinity, however strong, is an incidental rather than main function of these PFam54 proteins. The large number of PFam54 genes existing in any single B. burgdorferi genome may target different innate-immunity proteins of a single host species or the same immune protein of a variety of host species. Genetic variability of the PFam54 gene array suggests that universally present PFam54 lineages such as BBA64, BBA65, BBA66, and BBA73 may be better candidates for the development of broad-spectrum vaccines or drugs than strain-restricted lineages such as BbCRASP-1.